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_ , Abstract. We summarize recent results on the evolution of unpolarized parton 

^— V I densities and deep- inelastic structure functions in massless perturbative QCD. Due 

fSj ' to last year's extension of the integer- moment calculations of the three- loop splitting 

functions, the NNLO evolution of the parton distributions can now be performed 
reliably at momentum fractions x ^ 10~^, facilitating a considerably improved 
theoretical accuracy of their extraction from data on deep-inelastic scattering. The 
NNLO corrections are not dominated, at relevant values of x, by their leading small-x 
terms. At large x the splitting-function series converges very rapidly, hence, employing 
results on the three-loop coefficient functions, the structure functions can be analysed 
at N^LO for x >10~^. The resulting values for as do not significantly change beyond 
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Q^ i NNLO, their renormalization scale dependence reaches about ±1% at N'^LO. 
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1. Introduction 



Precise predictions for liard strong-interaction processes require transcending the 
standard next-to-leading order (NLO) approximation of perturbative QCD. Resumm- 
ations of large logarithms may be sufficient for specific processes, but generally full 
next-to-next-to-leading order (NNLO) calculations are called for. For electron-proton 
scattering and proton-(anti-)proton colliders this demands both partonic cross sections 
and parton distribution of NNLO accuracy. 

The former quantities are presently available only for the structure functions in 
deep-inelastic scattering (DIS) DH - which provide the backbone of our knowledge of the 
proton's parton densities and are among the quantities best suited for measuring Ug - 
and the total cross section for the Drell-Yan process - which in the form of W and Z 
production is an excellent candidate for an accurate luminosity monitor at Tevatron 
and the LHC 0. The calculation of other processes like jet production at NNLO is under 
way, e.g., the required two-loop two-to-two matrix elements have been computed [^ 
using the pioneering results [^ for the scalar double-box diagrams. See refs. for 
recent brief overviews. Partial NNLO results (the soft- and virtual-gluon contributions) 
have also been obtained for Higgs production via gluon-gluon fusion in the heavy top- 
quark limit 0, see also ref. 0. 

The three-loop splitting functions entering the NNLO evolution of the parton 
distributions have not been completed so far either |p. However, previous partial 
results 



10, 11 



have been substantially extended by the calculation of two more 
Mellin moments [|1^]. In Section 2 we discuss the resulting improvement |jl3[ of our 
approximations of the splitting functions in x-space WM , and compare, in the extended 
range x ^ 10~^ of safe applicability, the resulting approximate NNLO fiavour-singlet 



evolution and its scale stability to the NLO results. We also briefiy re-address [|I^, |r^ 
the question to what extent the leading small-x contributions to the splitting functions 
and coefficient functions dominate the small-x evolution. 

Taking into account the fast convergence of the splitting-function series shown 
in Section 2, the next-to-next-to-next-to-leading order (N^LO) corrections for the 
DIS structure functions can be effectively derived at x > 10"^ using available partial 



results |T^, |T2|, |18[ for the three-loop coefficient functions. The effect of the NNLO and 
N'^LO terms is discussed in Section 3 for the scaling violations of the non-singlet structure 
function F2 and the resulting determination of a 
predictions of the principle of minimal sensitivity 



191. Here we also illustrate the 



the effective charge method [^ 



and the Pade summation |22] which in this case, unlike the soft-gluon resummation, 
seem to facilitate a reliable estimate of the corrections even beyond N^LO. 

A first study has been performed |^ of the effects of the NNLO corrections (using 
our original approximations WM for the three-loop splitting functions) in a global parton 



analysis. See refs. |2^, ^ for beyond-NLO ag analyses of DIS data using methods more 
directly based on the integer- moment results |]TU|, 0. 



2. Singlet parton densities and structure functions at NNLO 

We first illustrate our approximation procedure for the three-loop splitting functions 
P^^'. As an example we discuss the N} term Pqg^i of the gluon-quark splitting function 
Pqg dominating the small-a: evolution of the quark densities. In the MS scheme employed 
in our studies, this function can be written as 

^9 i(^) = E ^- ln'"(l - ^) + /smooth (x) + 5: i?„ In" X + ^ Cp . (1) 

m=l n=l p=0 

The leading small-x coefficient Ci has been derived by Catani and Hautman [|1^. The 
function /smooth collects all contributions which are finite for < x < 1. This regular 
term constitutes the mathematically complicated part of Eq. (|l]), involving higher 
transcendental functions like the harmonic polylogarithms p6| . 



For our improved approximations [|I^ we choose three or two of the large-x 
logarithms in Eq. (|l|), a one- or two-parameter smooth function (low powers or simple 
polynomials of x) and two of the small-x terms {x~^ together with Inx or In^ x). Their 
coefficients are then determined from the six even-integer Mellin moments 

P!,a\iN) = fcixx^'-'PlfAx) (2) 



computed by Larin et al. [0 and Retey and Vermaseren [|l^]. By varying these choices 
we arrive at about 50 approximations (see Fig. 1 of ref. |l^). The two functions spanning 



the resulting error band for most of the x-range are finally selected as our best estimates 
for Pqg'i{x) and its residual uncertainty. 

These two functions, denoted by 'A' and 'B', are shown in Fig. 1 together with their 
(practically indistinguishable) real moments (|]) for 2 ^ A^ < 30 and their convolutions 

«|'«9lW=/'^/^?(!/)9(^) (3) 

with a typical gluon distribution g. In Figs. 1(b) and 1(c) the corresponding results 
for the Nj term have been supplemented for Nf = 4. Note that, like refs. |l^, |l^, we 
use the small expansion parameter as = a5/(47r); scaling down the ordinates by a factor 
2000 yields the results for an expansion in a^- The large impact of the A^ = 10 and 12 



moments [12| is illustrated by also showing our less accurate, but compatible original 



approximations ||I4| based on the four lowest even-integer moments . 



Knowing the leading x ^ In x term [|TT| is clearly instrumental in constraining the 



(2) 

small-x behaviour of Pqg^i - something not efficiently done by a small number of A^ > 2 
integer moments (|^). However, even at x < 10~^ where the non-x~^ parts contribute less 
than 10% to both approximations 'A' and 'B', this term does not sufficiently dominate 
over the (so far uncalculated) subleading CqX"^ contribution in Eq. (^, leaving us with 
a sizeable uncertainty of Pqg^i for x ^ 10~^. We will examine the dominance of the 
x~^ Inx and x~^ terms for the convolution (|^), which in any case considerably smoothes 
out the oscillating differences of the approximations, at the end of this section. 

After applying analogous procedures to the other three-loop splitting functions (see 
Figs. 2 and 3 of ref. |T^) we are ready to exemplify the effect of the NNLO contributions 
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Figure 1. (a) Exact results ||l^, Q (points) and approximations Q, Q (curves) 
for tlie moments of tlie N\ term of the tliree-loop gluon-quark splitting function in the 
MS scheme, (b) The approximations in x-space for N f = 4, and (c) their normalized 
convolutions with a typical gluon distribution. 



on the evolution of the parton densities. This is done in Fig. 2 via the logarithmic 
factorization-scale derivatives / = (iln//(iln/Xj of the typical singlet quark (S) and 
gluon [g) distributions 
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In Fig. 2(a) the resulting relative NNLO effects are shown for the standard choice 
[ir = fif of the renormalization scale. Note that the spikes close to x = 0.1 do not 
represent large NNLO corrections, instead they derive from zeros of the denominators. 
In fact, the NNLO corrections are very small (except for the quark evolution at small x), 
much smaller than their NLO counterparts not shown in the figure (except for the gluon 
evolution at small x, where already the NLO corrections are rather small). Consequently, 
as illustrated in Figs. 2(b) and 2(c), the renormalization-scale stability of the prediction 
is considerably improved over the full x-range. If measured over the conventional range 
1/2 Ai/ < fJr < 2/i/, the ^r dependence of E amounts to less than ±2% at large x and 
±5% small X, that of g to less than ±1% and ±2%, respectively. 

Also displayed in Fig. 2 are the present inaccuracies ('A' vs. 'B') of the NNLO 



results caused by the remaining uncertainties [T^\ of the three-loop splitting functions. 
These inaccuracies are entirely negligible at large x. Down to x ~ 10~^ they amount 
to about ±2% or less with respect to the central results not shown in the figure, even 
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Figure 2. (a) Size and present approximation uncertainties of the NNLO corrections 
to the scale derivatives of the singlet quark and gluon densities for the input (|^) and 
(pi) at iir — fJ-f- (b,c) The renormahzation scale dependence at NLO and NNLO for 
three typical values of x. 



if the bands in Fig. 2(a) are increased by 50% in order to account for any possible 
underestimate of the uncertainties. 

At lower scales the splitting-function uncertainties have a larger impact, mainly due 
to the larger value of a^. For example, the spread corresponding to Fig. 2(a) reaches 
about ±4% for S and ±3% for ^f at x = 10""^ and /ij ~ 3 GeV^ corresponding to as = 0.3. 
In view of the also enhanced NLO scale dependence, the approximate NNLO evolution 
represents an improvement over the NLO treatment even with inaccuracies of this size. 

The electromagnetic singlet structure function F2 and its Q^ derivatives are 



presented in Fig. 3 for the parton densities (HI) at Q^ 



^/,o 



30 GeV^ The large 



NNLO corrections at very large x originate in the non-singlet part of the two-loop 
quark coefficient function. Note, however, that the (positive) gluon contribution to 
dF2,s/d\nQ'^ at NNLO still amounts to 5% at x = 0.5 (40% more than at NLO), and 
falls below 1% only above x = 0.7 |]T3[. This effect is large enough to jeopardize analyses 
applying a non-singlet formalism to the proton's F2 in the region x > 0.3. 

The negative NNLO corrections to F2 at small x arise from the two-loop gluon 
coefficient function C2J. The Q^ derivative, on the other hand, receives a +10% NNLO 
correction for 10~^ ^ x ^ 10~^ ; its break-up is illustrated in Fig. 3(a) by the results for 
P*^^^ = 0. Also for F2 and its scaling violations the inclusion of the NNLO terms leads to 
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Figure 3. (a) The NNLO corrections for the singlet F2 and its Q^ derivatives (linear 
at small x, logarithmic at large x) for the input (^ at fif — ^y = fi^ = ^"^^-^ = Q^ « 
30 GcV^. (b,c) The scale dependence at NLO and NNLO for three typical values of x. 



obtained by keeping only the x ^ In a; terms 
with the (for P*^^) approximate) full results. The 



a substantial decrease of the scale uncertainties as shown in Figs. 3(b,c), which facilitates 
more precise extractions of the parton distributions from data on these quantities. 

We conclude this section by examining the dominance of the small-x terms of the 
NNLO splitting functions and coefficient functions for the small-x convolutions. In Fig. 4 
the results for f®g, f = P^J>, P^f i 
of / are compared, down to x = 10 
dependence on the gluon distribution g is illustrated by employing, besides our 'steep' 
standard input (|), also a low-scale 'flat' shape, xg ~ x° for a; — i> 0. 

Keeping only the leading x~^ In x terms |]TT| of P^j^ and Pj^^ does not lead to 
reasonable approximations as shown in Figs. 4(a,b), regardless of the gluon distribution. 
Even for the more favourable fiat shape the offsets amount to about a a factor of two 
even at x = 10~^. Besides the x~^ contributions, the non-x~^ terms do not seem to be 
sufficiently suppressed either, at least for a steep gluon distribution. Due to the present 
large uncertainties [0 on the x~^ terms, however, definite conclusions especially for a 
flat gluon distribution require the computation of the exact x-dependence of P/j . 
On the other hand, such conclusions can be drawn already |]15| for the convolutions of 



the two- loop coefficient function C2J [IH shown in Fig. 4(c). The leading x ^ term [^ 



does not dominate over the non-x ^ contributions at any x- values of practical interest. 
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Figure 4. The convolutions of the leading (and, for p(^\ subleading) small- 

X terms of (a,b) the NNLO splitting functions Pf^ and Pf^ and (c) the NNLO 
coefficient function c\ with two typical 'steep' (upper row) and 'flat' (lower row) 
gluon distributions, compared to the respective (for P*^^' approximate) full results. 



3. Non-singlet structure functions at NNLO and beyond 

The scaling violations of the non-singlet structure functions -Fa,NS, ^ = 
conveniently discussed in terms of the physical evolution kernels -ft'a,NS ■ 
d „ 



a,NS 
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cilnQ2 
The N'LO expansion coefficients Ka^i are combinations of the coefficient functions up to 



/ loops and the splitting functions up to Z+l loops. See Eqs. (2.7) - (2.9) of ref. [|T9[ for 
the details. The advantage of Eq. (|^) is that any dependence on the factorization scheme 
and the scale /xj has been eliminated explicitly. Note, however, that in a^ analyses in 
terms of coefficient functions and parton distributions the main uncertainty arises from 
the dependence on the renormalization scale /x^, not the factorization scale /i/. 

The NNLO kernels Ka^2 are fixed by the two-loop coefficient functions |l[] (for which 
compact parametrizations are available |]n[) together with our approximations [^ for 
three-loop splitting functions. The uncertainties of the latter are negligible in the region 
X > 10^2 in which we are mainly interested for the non-singlet case. 

Due to the fast large- a; convergence of the splitting function series illustrated in 
Fig. 2, the N^LO kernels are dominated by the coefficient functions. The seven lowest 
even (for -^1,2) and odd (for F3) moments at three loops have been computed by Larin 
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Figure 5. The PIMS, ECH and Fade estimates of (a) the NNLO, (b) the N^LO, and 
(c) the N^LO contributions to the A^-space evolution Icernel for iF2_NS a-t Mr = Q^^ 
Nf = 4: and as = 0.2. Also shown are the (approximate) full NNLO and N'^LO results. 



et al. |[T^ (A^ < 10 for -^1,2) and Retey and Vermaseren [Q. When complemented by the 
four leading large-x terms (1 — x)~^ In {1 — x), k = 2, . . . , 5 derived |T^ from the soft- 
gluon resummation |]T^ , this information facilitates x-space approximations (analogous 
to those for P^^-* exemplified in Section 2) which are sufficiently accurate in the above- 
mentioned region of x, as shown in Figs. 2 and 3 of ref. [|l^]. The small effect of the 
uncalculated four-loop splitting functions P^^^ can be estimated by assigning a 100% 



,(3), 



,(3)/ 



error to the Fade prediction P^s (^) ~ -^ns (^) [1/1] Fade • Actually the small residual 
uncertainties of Ka^^ at x > 10~^ are dominated by this error, not by the approximation 
uncertainties [^ of the three-loop coefficient functions. 

Knowing the evolution kernels to such a high order also facilitates a test of the 
predictions for Ka^i derived from the Fade summation ^2\ of the perturbation series 
and from renormalization-scheme optimizations like the principle of minimal sensitivity 
(FMS) pg and the effective charge (ECH) method ^. In Figs. 5(a) and 5(b) these 
predictions are compared to the (approximate) full kernels i^2,2 and 7^2,3 for -F2,ns in 
A^-space. The estimates by all these methods agree rather well for the a^ (N^LO) 
contribution shown in Fig. 5(c). Such an agreement is usually interpreted |^ as 
evidence for the approximate correctness of the predictions. 

The large- A^ / large-x behaviour of the kernels -ft'a.NS is dominated by the soft-gluon 
terms a^,"^^ In N, 1 < k < l + l. The soft-gluon resummation |T^ at next-to-next-to- 
leading logarithmic accuracy [0 fixes the coefficients of the leading three terms at all 



orders. However, as shown in Figs. 8 and 9 of ref. |T^, these contributions do not provide 
quantitative predictions for Ka^i at practically relevant values oi N / x, due to the (quite 
generally, see ref. [^) large coefficients of subleading logarithms at higher orders [|il 
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Figure 6. (a) The pcrturbative expansion of the logarithmic Q^ derivative of F2,ns 
for the input (0) for Nj ^ A and ^l = Q"^ . (b) The N"LO contributions, n = 1 . . . 4, 
to the results shown in (a), and (c) their renormalization scale uncertainties. 



The effect of the higher-order terms (using the Mellin inverse of Fig. 5(c) at N^LO, 
see Fig. 12 of ref. |T3 for the Fade estimates of the corrections beyond N^LO) is 
exemphfied in Fig. 6 for the logarithmic Q^ derivative of 



F2Ns(a;,Qo - 30 GeV^) = x'^-\l - xf , as{Ql) = 0.2 . (7) 

Also shown are the renormalization-scale uncertainties estimated using the conventional 
interval \Q'^ < fi^ — 4Q^. The a^ (N^LO) corrections and the N^LO scale dependence 
are very small at a; < 0.6 ; they both reach about 2% and 5% of the total results only at 
X = 0.65 and 0.85, respectively. The corresponding numbers at N^LO read 1% and 3%. 
Sample fits of a^ to the Q^ derivatives at 0.05 < x < 0.75 for Q^ = Ql (using the 
N^LO predictions of Eq. (|^ as model data) yield |TD 



as(<5, 



o)nlo 



0.2080 



+ 0.021 

- 0.013 

+ 0.003 

- 0.001 



'^s(vo)nnlo 



0.2010 



+ 0.008 
- 0.0025 



«.(Q^)n3lo = 0.2000 ^ 121 , aMULO = 0-2000 ^ ^^o, ' ^^^ 
where the errors include the above /i^ variation and, at N^LO and N^LO, the small 



approximation uncertainties. Corresponding results for F^ can be found in ref. [19 
In both cases the N^LO and N^LO corrections, unlike the NNLO terms, do not cause 
significant shifts of the central values, but just lead to a reduction of the fir uncertainties 
which reach about ±1% at N^LO. 

Fortran subroutines of our approximations of splitting functions |T^ and non- 
singlet coefficient functions [IS] at three loops, and of the parametrizations [Q of the 
two-loop coefficient functions [0 and the convolutions entering the physical evolution 
kernels |T^ can be found at pittp : //www . lorentz . leidenuniv . nl/| ^ avogt. 



10 

Acknowledgment 

This work has been supported by the European Community TMR research network 
'QCD and particle structure' under contract No. FMRX-CT98-0194. 

References 

n 



[2] 
[3] 

[4; 

[5] 
[6] 

[7 

[8] 
[9] 

[lo; 
[11 



[12; 
[13; 

[14 

[15; 
[16; 



[1?: 



[is; 

[19 

[20' 
[21 
[22; 

[23; 
[24 

[25; 
[26; 
[27; 



E.B. Zijlstra and W.L. van Neerven, Phys. Lett. B272 (1991) 127; ibid. B273 (1991) 476; ibid. 

B297 (1992) 377; Nucl. Phys. B383 (1992) 525 
R. Hamberg, W.L. van Neerven and T. Matsuura, Nucl. Phys. B359 (1991) 343; E.B. Zijlstra and 

W.L. van Neerven, Nucl. Phys. B382 (1992) 11 
M. Dittmar, F. Pauss and D. Zurcher, Phys. Rev. D56 (1997) 7284; V.A. Khoze, A.D. Martin, 

R. Orava and M.G. Ryskin, Eur. Phys. J. C19 (2001) 313; W.T. Giele and S.A. Keller, 



hep-ph/0104055 



C. Anastasiou, E.W.N. Glover, C. Oleari and M.E. Tejeda-Yeonians, Nucl. Phys. B601 (2001) 



318; ibid. B601 (2001) 341; |hep-ph/0101304| ; |hep-ph/0104178 



V.A. Smirnov, Phys. Lett. B460 (1999) 397; J.B. Tausk, Phys. Lett. B469 (1999) 225 



W.L. van Neerven, Pramana 55 (2000) 101 (|hep-ph/0003043| ); E.W.N. Glover, |hep-ph/0106069 



S. Catani, D. de Florian and M. Grazzini, JHEP 05 (2001) 025; R. V. Harlander and W. B. Kilgore, 

Phys. Rev. D64 (2001) 013015 
M. Kramer, E. Laenen and M. Spira, Nucl. Phys. B511 (1998) 523 
S. Moch and J.A.M. Vermaseren, Nucl. Phys. (Proc. Suppl.) 89 (200) 131, 137; S. Moch, J.A.M. 

Vermaseren and M. Zhou, work in progress 
S.A. Larin, T. van Ritbergen, and J.A.M. Vermaseren, Nucl. Phys. B427 (1994) 41; S.A. Larin, 

P. Nogueira, T. van Ritbergen, and J.A.M. Vermaseren, Nucl. Phys. B492 (1997) 338 
S. Catani and F. Hautmann, Nucl. Phys. B427 (1994) 475; J. Bliimlcin and A. Vogt, Phys. Lett. 

B370 1996) 149; V.S. Fadin and L.N. Lipatov, Phys. Lett. B429 (1998) 127; M. Ciafaloni and 

G. Camici, Phys. Lett. B430 (1998) 349 
A. Retey and J.A.M. Vermaseren, Nucl.Phys. B604 (2001) 281 
W.L. van Neerven and A. Vogt, Phys. Lett. B490 (2000) 111 
W.L. van Neerven and A. Vogt, Nucl. Phys. B568 (2000) 263, ibid. B588 (2000) 345; see also 

Nucl. Phys. (Proc. Suppl.) 89 (2000) 143 
W.L. van Neerven, talk at the 1993 DESY Workshop (unpubhshcd) 
J. Bliimlcin and A. Vogt, Phys. Rev. D57 (1998) 1; ibid. D58 (1998) 014020; J. Bliimlcin, 

V. Ravindran, W.L. van Neerven and A. Vogt, Proceedings of DIS 98, Brussels, April 1998, 



eds. Gh. Coremans and R. Roosen (World Scientific 1998), p. 211 (liep-ph/9806368) 
G. Sterman, Nucl. Phys. B281 (1987) 310; L. Magnea, Nucl. Phys. B349 (1991) 703; S. Catani and 

L. Trentadue, Nucl. Phys. B327 (1989) 323; ibid. B353 (1991) 183; S. Catani, G. Marchesini 

and B.R. Webber, Nucl. Phys. B349 (1991) 635 
A. Vogt, Phys. Lett. B471 (1999) 97; ibid. B497 (2001) 228 
W.L. van Neerven and A. Vogt, Nucl. Phys. B603 (2001) 42 
P. M. Stevenson, Phys. Rev. D23 (1981) 2916 
G. Grunberg, Phys. Rev. D29 (1984) 2315 
M.A. Samuel, J. ElHs and M. Karhner, Phys. Rev. Lett. 74 (1995) 4380; J. Ellis, E. Gardi, 

M. Karliner and M.A. Samuel, Phys. Lett. B366 (1996) 268; Phys. Rev. D54 (1996) 6986 
A.D. Martin, R.G. Robert, W.J. Stirling and R. Thorne, Eur. Phys. J. C18 (2 000) 117-126 
A.L. Kataev, G. Parente and A.V. Sidorov, Nucl. Phys. B573 (2000) 405; |hep-ph/0012014 ; 



hep-ph/0106221 



J. Santiago and F.J. Yndurain, Nucl. Phys. B563 (1999) 45; tiep-ph/0102312 



E. Remiddi and J.A.M. Vermaseren, Int. J. Mod. Phys. A15 (2000) 725 



E. Gardi and J. Rathsman, liep-ph/0103217 



